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Abstract 

The importance of studying matter at high p increases as more 
astrophysical data becomes available from recently launched space- 
crafts. The importance of high T studies derives from heavy ion data. 
In this paper we set up a formalism to study the nucleons and isobars 
with long and short range potentials non-pertubatively, bosonizing 
and expanding semi-classically the Feyman integrals up to one loop. 
We address the low density, finite T problem first, the case relevant 
to heavy ion collisions, hoping to adresss the high density case later. 
Interactions change the nucleon and isobar numbers at different p and 
T non-trivially. 
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1 Introduction 



Isobars (A's) play a very important role in nuclear physics. Two pion ex- 
change with a A intermediate state is known to produce the intermediate 
range attraction between nucleons[p. Many attempts have been made to 
extract the effect of A intermediate states in nuclear matter dynamics, and 
microscopical calculations have shown the overwhelming relevance they have 
in binding of nuclear matter ||l|, 

A is the next excited state in the nucleon spectroscopy with a large de- 
generacy factor (4x4 for spin and isospin). Paradoxically however not too 
much efforts have been devoted to extract the amount of A population in the 
nuclear matter ground state 0, It is on the other hand conceivable that 
the small numbers found so far (Ref. ^ clearly overestimates the A compo- 
nent of the nuclear ground state) suggest that the presence of A's could at 
most induce rather small effects in the ground state properties. Only recently 
use of modern A^-A^ interaction has renewed the interest in this field [Q, but 
still dynamical aspects like spectral functions and response functions deserve 
further investigations even taking advantage of the modern computing facil- 
ities. 

Even disregarding, however, the relevance of traditional nuclear matter 
calculations, new achievements both on the astrophysical sector and on the 
side of ultra-relativistic heavy-ions collisions suggest to explore the properties 
of the hadron matter at finite (high) temperature. One should keep in mind 
in fact that in the latter low density hadronic matter is expected to 

be the relic of a phase transition to and from a quark- gluon plasma phase. 
The knowledge of the nuclear matter properties at low density and high 
temperature is thus linked (maybe not in a simple way) to the signals that 
such a phase transition has indeed occurred. But such conditions seem to be 
much more sensitive to the presence of A's, as real A's can be produced. 

This has already been established by calculations where non-interacting 
nucleon and isobar excited states are embedded in a thermal bath ( see for 
example Bebie et al. and Dey et al. [^). Interactions between nucleons 
(A^ in short) and A's being strong, it is very difficult to handle. A non- 
perturbative treatment is mandatory. 

The framework we use in this paper to deal with a system of interacting 
nucleons and isobars is provided by the path integral technique developed by 
one of us and his co-workers p, ^ |lT|. It involves Bosonisation arising 
from integration over the fermionic fields. 
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For non-interacting fermions the methodology is simple. However, for 
long range one pion exchange potential ( OPEP ) or short range correlations 
(SRC in short) non-perturbative semi-classical approximation schemes are 
necessary. This provides a mean field. Higher order corrections can be put 
onto it. 

The method is described in short in the next section. The nucleon and 

the isobar, being heavy particles can be treated non-relativistically. This 
allows the Fermi integrations and the Pauli principle to be dealt with. As 
the density , p, increases the baryons are packed more closely so that the 
Pauli principle is expected to play a crucial role. 



2 The partition function for the Nucleon in- 
teracting with the Pion 

To exemplify, we consider a system of nucleons interacting with pions. OPEP 
potential between nucleons and pions is written as 

v;(q) = rf )v.(q)rf ) (1) 

where 

V.(q) = - ^'T 2^ 2 (2) 

and 

rf ) = <Ti • qr, . (3) 

The partition function in the grand canonical ensemble, in terms of path 
integrals, reads 

where H , /j, and are the Hamiltonian , the chemical potential and the 
number operator respectively. 

For non-interacting nucleons we can easily evaluate the partition function 
in terms of the Fourier transforms for the fields ^'(q, ^n) '■ 

Zo = Z = j D[^p\ ^]eE„/<i'9V't{-q,-c^n)6o-'{9>'^")V'{q,a;n) (5) 
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where QQ{q,Un) is the free Green's function at finite temperature : 

Qo{q,^n) = 



(6) 



2M 



fl - lUJn 



with 



(2n+ l)7r 



(7) 



the functional integration becomes the determinant of the Green's function. 

To evaluate the partition function with OPEP interaction , eqn.(^), an 
approximation scheme of boson loop expansion is used. This scheme has been 



explained in details in the already quoted refs. P, |TT|. Basically, func- 
tional integral over the fermion fields leaves an effective bosonic Lagrangian 
which in the zeroth order gives rise to a mean field like that of RPA and in 
the 1st order to the one boson loop approximation. 

This is achieved by availing of the Hubbard-Stratonovitch transformation 



12| , pT^ , |Tj] to convert the potential part inside the Feynman integrals into 
a functional integral. For example, the potential inside eqn.(^), under the 
Hubbard-Stratonovitch transformation, becomes 



D[(j]e^ /o^f^T / d? X d'^ytJi{x)V:;;^ {\x-y\)a,{y)+i dr f d^xai{x)tp^x)ritp{x) 



d(Ji{xk) 



where the measure over a,- is defined as 



Performing the functional integration over the fermionic fields it becomes 



Z = Vdet(-V.) j D[(Ti]e2-'^--» det [G^^ - aiPj . (9) 
Since Qq^ and cTiFj commute : 

Z = ^det(-V^)det j /}[o-,]e^'^>^"'^'+T'^i°s[^-^o'^'^< (10) 
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3 Introduction of A's, p-mesons and Short 
Range Correlations (SRC) 

The simplest way to introduce the A's in this formahsm is to generahse the 
nucleon field to a column vector 




and the Fj (eqn.|D as 



cr • qTi 



I /zEiVA gt . qTt Uaa^ . 7: 

\ JttNN JttNN 



■ qT, 

JttNN 
/ttZ 



[12) 



Free Q% becomes 



g%{q,UJn) 



1 



+ 6M — fiA — i^n 



(13) 



2Ma 



where 6M = Ma — M. Matrices Si and % are the generators of the 3/2 
representation of the SU{2) group ( see for instance [|lll). But Si and Tj are 
the familiar 4x2 transition spin - isospin operators (see for example, ||15|| ). 

Further, a reasonable description of the dynamics of a nuclear system 
interacting via meson exchange needs to account, in order to be realistic, of 
short range correlations (SRC). This SRC has been introduced in nuclear 
systems in many ways, for example through the Landau-Migdal parameter 
It has been shown by 
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that at the one boson loop level, beyond 
the Landau-Migdal theory, a momentum dependence of g' is necessary. Now 
there are two independent parameters, namely the value of the function g'{q) 
at the origin and a cut-off that kills the effective interaction due to nucleon- 
nucleon (or nucleon-A) short range repulsion. Both parameters are density 
dependent, in a way which is not known. 

A simpler treatment to deal with SRC is given by Brown et al 



which is successfully applied by Oset and coworkers (see, among many other 
papers, [|17|)- It amounts to multiply the nucleon-nucleon interaction V(r) 
by the two-body correlation function g{r), i.e.. 



V(r) — > g{r)V{r) 



(14) 
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and further to approximate g{r) with 



air) ~ -joiqcr) (15) 

where qc is taken roughly equal to the tu-mass. This approach not only works 
well but also is almost parameter free. Density dependence of the correlations 
is automatically taken care of. It is to be noted that, it amounts to fix g' at 
about 0.7. 

Of course a A-hole pair can be excited by a p-meson as well. We shall see 
in the following that a p-meson exchange, with the same scheme for SRC as 
before, can also be included in the formalism. In the following calculations, 

= 770 MeV/c and qp = 1200 MeV/c. 



4 The Loop Expansion 



We now expand the partition function (eq.|T0D semi-classically. At the saddle 
point, i.e., at the mean field level, Uj = and the partition function becomes 



^mean field 



det — Vtt det Qq ^ . 



(16) 



Quadratic deviations of the field from its mean field value lead to 

L>[(T,]e^'''^-"''''-^^'""^'^'"'^^^^« = {det V^^ - Trri^orj^o}^' (17) 

so that the partition function (eq. |10|) at the one loop order becomes 

Z = det ^-ie'^''^°^[^"'^'^^'^°^^^°] = Zoe^'^°^^^~^'''^'^''^'^°^ (18) 

where 



n 1+ 



2AfA 



+5M~flA 



(19) 



The grand potential, Q = logZ, in one-loop order becomes 
= ~ {logZo + Tr log [/ - V.T,,goT,go]} 



(20) 
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Thus for the 0-order part we get 



^0^ — 



V f d^q 



P J (27r)= 



log 



1 + e 



+ log 



1 + e 



+SM-fiA 



(21) 

Derivatives of Q with respect to /in and //a result in N and A numbers 
respectively. Therefore, 



V 



where 



<Na> = V 
n%{q) = : 



d^q 
(2^ 



nliq) 



1 + e 



1 



1 + e 

For first order correction we define 

1 







(22) 
(23) 

(24) 
(25) 



(26) 



(the trace makes i j vanishing) which turns out to be the generalisation of 
the Lindhard function to the finite temperature case. We have three different 

n° — s, corresponding to the N particle - N hole, the A particle - N hole 
and the A particle -A hole propagators. The last one vanishes at the T ^ 
limit. 



n^A(q,i^n) 



+ 



I 



d^k n%{k) - n%{q + k) 



(27) 



(27r)3 + e(fc) -e{q + k) 

16 r _fk n%{k) - n%{q + k) 

9 J (27r)3 iun + e{k) - e^{q + k) + /ia - fJ'N 
{N ^ A) (28) 



16 



d^k <(A;) - nl{q + k) 
(27r)^ iun + eA{k) -eA(q + k) 



(29) 
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where the numbers in front are contributions from the spin-isospin traces and 



/9 



(30) 



The parity of these function with respect to simphfy them further. 
Thus the total 11 is given by 

n = f^NN^%N + fnNA^%A + fnAA^^AA ■ 

SO that one loop grand potential (eq.^OD is 



(31) 



f]i-i°op = Tr log [i + v^q'n] =no + -J2-TT [v^nr (32) 

P P n=l ^ 

1 

Using the identity ^ = / X"'~^dX and substituting = XVj, we get 



n 



1— loop 



Further, summation of the series results in 



Tr 



(33) 



(34) 



The remaining trace means integration over the 3-momentum and sum over 
the frequencies. 

Number of A's and nucleons can now be obtained easily. 



< Na >=< Na >° -^Tr 



/3 " 1 - v^n di^A 



and 



< A^^ >=< A^^ >o -^Tr 



1 ^ K dU 



(35) 



(36) 



13 1 - V^n d^iN 

The 0-order quantities are given by (^) and (p3D . To evaluate eqs. (p5| 
and (BBf) we find that Ist-order correction term (n=l of eq. (P^D) 



1„ d 



< Nn >^''-°''^= -^Tr 



[v.n] 



needs added attention because it is ill defined as in the case of self-energy of 
the nucleon in a medium : 

S(A;,c..) = ^E / ^^o(q,^n'R(|k-q|). (37) 

We need a convergence factor e*'^""' which produces an extra factor n^^{q). 
The derivative with respect to fi^ gives 



< >lst-ord^ 



■IN 

12 



(38) 



^'S.(lk-ql) ^ 



16 



(27r)3 
d^k d 



n%{q)n%{k) 



(27r)3y (27r)3 9/i^ 
Similarly one gets for the A contribution 



n%{q)n^{k) + n%{k)nA{q)] V,(|k - q|) 



16 

y 



>= 

d^q 



(39) 



d^k d 



(27r)3y {27T)^dfiN 

d^q f d^k 



n%{q)nA{k) + n%{k)nA{q)] V^(|k - q|) 

The remaining part of the series can be evaluated straightforwardly : 

d^q v^n dU 



(27r)3l-V^n9/x^,A ■ 



(40) 



by evaluating 11° -s (eqs. (p^^-pop. 

Some remarks are not out of place here. 



1. Sum in eqn. ( ^OD can be further transformed into an integral over the 
imaginary axis. But this does not make numerical calculation easier. 

2. the present formalism amounts to evaluate the diagram of fig. |I] where 
the black dot denotes the dU/dfi insertion and the dashed circle the 
RPA-dressed pion. Since pions are spin-longitudinal while the p's are 
spin-transverse and in an infinite medium 11 conserves the elicity, the 



9 



two kinds of mesons cannot mix together. Thus the p meson contri- 
bution simply amounts to the similar diagram with the pion exchange 
potential replaced by the p exchange with a factor 2 in front of 11 
coming from spin traces. 

3. In general, |^ in eqn. (|40D obtained from eqns. (p6|-^) becomes ill- 
defined for n = 0. But if eqns. (p6| - ^0D describe analytical functions, 
even for n = 0, the singularity can be separated from the integral and 
handled analytically. 



5 Numerical results: densities and chemical 
potential 

The previous section provided the tool needed for the evaluation of < > 
and < A^A >• At the equilibrium Pn = Pa = P and hence both < > 
and < A^A > can be regarded as functions of p. Putting for short < >= 
Nj^{p) and < A^a >= N/^{p) it is also clear that 

Nn{p) + Na{p) =A = Vp = V— I , (41) 

where A is the baryon number and the last equality can be regarded as a def- 
inition of an effective Fermi momentum. Here p is the input of the problem. 
If instead we choose A in order to fix the thermodynamical conditions of the 



system we need to consider (^) as an equation for p and solve it numerically. 

Fig. 1^ displays the results of such a calculation. It shows the behaviour 
of the chemical potential at the equilibrium for different values of the tem- 
perature and the nuclear density. The graph contains only pion exchange at 
three different levels of complexity, namely calculations up to second order 
(dashed lines), full one- loop calculation (solid line) and full calculation with 
dynamical pion (dotted line). We remember that in the frame of Matsubara 
Green's function the pion propagator reads 

1 

A, 



The figure displays five bunches of lines, corresponding to five different den- 
sities, namely (starting from the above), p/po = 0.2, 0.4, 0.6, 0.8, 1; of course 
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Po denotes the normal nuclear density. Remarkably these three different dy- 
namics are almost completely superimposed: only looking very carefully at 
the figure one sees that the dashed line is a little bit higher than the other 
two, that instead are practically coincident. 

The second plot we present (fig. ^ displays instead the ratio between 
nucleon and A density keeping fixed the sum Pn + Pa = P again as a function 
of the temperature and with five bunches as before. Of course this time the 
highest bunch corresponds to the highest density. Here the plots are much 
more sensitive to the dynamics. 

Insofar we have only considered exchange of pions. Let us add some 
more dynamics by including the exchange of correlated p-mesons the results 
are shown in figs. ^ and ^. Remarkably, the p-meson exchange lowers 
significantly the chemical potential, while the ratio pa/pn as a function of 
the temperature is sensible to it only in the low T region, where the full 
dynamical calculation shows a remarkable reduction of the pa/pn ratio. At 
higher T instead the effect fades out. 

The model at hand contains only a few number of parameters. Coupling 
constants and cut-off are more or less well defined (we use here standard 
coupling Cp = 2.3) the only crucial parameter being the one connected to 
SRC, and in particular q^. Thus in figs. |] and ^ we plot again the temper- 
ature and the ratio Pa/pn for the most complete dynamical case (vr plus p 
fully dynamic at the one- loop order but with different values of g^, namely 

= 650 MeV/c, (solid line), q^ = 770 MeV/c (dashed line) and q^ = 900 
MeV/c (dotted line); the same conventions as in the two previous figures are 
adopted. We observe that while the chemical potential is almost insensitive 
to the variations of the ratio pa/pn, on the contrary, shows that this ra- 
tio sensibly increases with g^, in particular in an intermediate range of the 
temperature, while at temperatures of the order of 150 MeV the effect seems 
to fade out. 



11 



6 Numerical results: the momentum distri- 
bution 



The momentum distribution displays an impressive dependence upon the 
temperature. We define the momentum distribution through the expression 

<N>^ = Trj-0^^n^{q) (42) 

< >A = Tr I -0^nM (43) 

the trace being taken over spin-isospin degrees of freedom (4 for nucleons, 16 
for A's). In practice this amounts to take the integrand of all the previous 
calculations. Thus the result is a byproduct of the previous calculation. The 
results are plotted in fig. ^ for different temperatures and densities. The 
normalisation is everywhere to a 6{kp — q) when T ^ 0. 

It is impressive to remark that the distribution is immediately spreaded 
to high momenta with the increasing of the temperature. Remind that the 
distribution is weighted with a in the Jacobian, so that the high momentum 
component become more and more important. An even more striking insight 
can be obtained by plotting the quantity q^n{q) and by comparing it with 
the Free Fermi Gas result at temperature, as we do in fig ^. There it is 
clear that the areas of the curves is preserved while the shape is enormously 
modified. 

Furthermore, at high temperature the dependence upon the density seems 
to become weak and something like an universal behaviour could be sug- 
gested. 



7 Discussions 

People speculate about the number of isobars in a hot soup of hadrons pro- 
duced in a heavy ion collision. Even assuming that the short time available in 
heavy ion collisions is enough to thermalize the baryons and produce a uni- 
form collective fiow, the dynamics discussed in this paper is not sufficient to 
enlight the great variety of processes that could occur at finite temperature. 

The vacuum at finite temperature is known to contain pions, since these 
particles are light. Dey et al |T^ have shown that finite temperature will 
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couple channels with different parity and isotopic spin. For example the p 
and the Ai meson along with a longitudinal pion mix to order T^. The poles 
do not move till the next order T'*. In the same way the nucleon at finite 
T will couple to the A1/2 state and the isobar to the nucleon excited state, 
N*j^{lb2Q) as well as the odd parity isobar A(1700) H]. 

Further, if we focus on the relic of a previously realized quark-gluon conde- 
sate, then the strangeness contribution should be relevant, and investigations 
on strange hyperons is also suggested . 



Coming to more conventional degrees of freedom, it was shown by Be- 
daque in a long letter, using chiral perturbation theory, that the nucleon 
mass increases a little, but that the nucleon acquires a substantial width and 
the mass is decreased so that the isobar-nucleon splitting becomes smaller. 
In chiral limit Bedaque's result would appear in order T^. 

The present considerations seem thus suggest as future perspectives the 
extension of the present 1-loop calculation by one side to the study of other 
observables in a hadron gas, like for instance effective masses and widths 
and also the entropy of the system; and on the other side to the extension 
of the present scheme to a richer dynamics encompassing vector mesons and 
strange hadrons. 

A further interesting development, to be carried out in the future, will be 
to extract, from the present formalism, the corresponding behaviour as an 
expansion in powers of the temperature. 
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Figure 1: One-loop diagram for the number of A and nucleons. 
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Figure 3: Ratio p\/ Pn as a function of the temperature. Lines as above. 
The total densities are 0.2 poi 0.4 Pq, 0.6 0.8 po and 1 po in ascending 
order. Only pions are accounted for. 
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Figure 5: Ratio Pa/ Pn as a function of the temperature, with pion plus p 
exchange. Lines as in ^ 
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Figure 7: Ratio Pa/ Pn as a function of the temperature, with pion plus p 
exchange. Lines as in |^ 
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Figure 8: Momentum distribution of nucleons (solid line) and A's (dashed 
line) at different densities and temperatures 
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Figure 9: Momentum distribution multiplied by (f of nucleons (solid line) 
and A's (dashed line) at different densities and temperatures. The dotted 
line corresponds to the Free Fermi Gas at temperature. 
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